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Method and apparatus for producing contrast agents for magnetic resonance imaging 

Field of the Invention 

5 The present invention relates to an apparatus and method for para-hydrogen induced 

hyperpolarization of a compound, and in particular for the preparation of a contrast agent for 
magnetic resonance imaging procedures. 

Background of the Invention 

Magnetic Resonance Imaging (MRI) is an important diagnostic technique. It is especially 
10 attractive since it is non-invasive and does not expose the patient to potentially harmful 
radiation such as X-rays or radiation from radioactive materials. Significant progress has 
recently been made in the quality of images as well as finding new applications of the 
technique. The progress relies on a rapid development of the digital image processing, the 
refined Nuclear Magnetic Resonance (NMR) techniques and the development of effective 
15 contrast agents (imaging agents). An emerging technique of particular interest involves 

Magnetic Resonance (MR) contrast agents based on the principle of pre-polarization of the 
nuclear spins, also called hyperpolarization 

For the resonance phenomena, which is the basis of NMR and MRI, to occur, isotopes with 
non-zero nuclear spin have to be present. In addition, since NMR is not an extremely sensitive 

20 technique, a relatively high concentration and/or a high gyromagnetic ratio is needed, 

especially for imaging purposes. The use of a selected isotope in a contrast agent which is to 
be injected in a patient, puts further requirements on the selected isotope, for example 
regarding toxicity. A number-of isotopes have the required spin properties, but only a few are 
considered interesting for the use in contrast agents, for example the carbon isotope 13 C and 

25 the nitrogen isotope l5 N. The carbon isotope 13 C has many properties that would it make it 
useful as a functional part of an MRI contrast agent. An important feature is the long 
longitudinal relaxation time, Tj. The relaxation time needs to be long in order to have time, 
after the generation of the contrast agent, to inject the contrast agent into the patient and to 
allow the contrast agent to be transported to the organ that is to be studied. To make a useful 

30 MR-contrast agent of this kind, the signal strength has to be boosted significantly over the 
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thermal equilibrium signal. Patent application WO 00/71 166, by the same applicant, describes 
a process and an apparatus for increasing the polarization and hence the signal from a small 
organic molecule containing for example 13 C. The signal was increased with a factor 10 4 . The 
process is referred to as Para-Hydrogen Induced Polarisation (PHIP) and may comprise the 
5 transferring of nuclear spin-order from para-hydrogen to spin polarization in non-zero spin 
nuclei in molecules, e.g. to ,3 C or ,5 N nuclei. 

Hydrogen molecules exist in four different spin states. In one of the forms, characterised by 
antiparallel spins, the magnetic moments of the protons cancel. This form is called para- 
hydrogen. The other three forms, with a net magnetic moment, are referred to as ortho- 

10 hydrogen. The para-hydrogen will not rotate at low temperature whereas the ortho-form must 
rotate with a high frequency at all temperatures because of quantum-mechanical symmetry 
requirements of the wave function. This indicates that at low temperature the para-form will 
have a significantly lower energy, and hence is the energetically favoured form. At 
temperatures below 20 K the equilibrium ratio of para- and ortho-hydrogen approaches 100:0, 

15 at 80 K the ratio is 48:52 and at room temperature approximately 1:3. The equilibration can 
be speeded up by the presence of a transition metal catalyst, e.g. Fe2Q3. Para-hydrogen relaxes 
slowly (if no catalyst is present) at room temperature. 

In WO 00/71 166 it is described how to catalytically hydrogenate (with para-hydrogen) 
unsaturated compounds comprising non-zero spin nuclei such as 13 C. The spin correlation of 

20 the protons from the para-hydrogen will be preserved during and after hydrogenation, and the 
influence on the spins of the l3 C nuclei breaks the symmetry of the spin system. The protons 
will now give an NMR-signal, but the non-equilibrium spin order is not sufficient to make the 
molecule useful for imaging purposes since it has an anti-phase behaviour that is not ideal for 
imaging. In the above cited applications, and further in 'Tarahydrogen-Induced Polarization 

25 in Imaging" by K. Golman et aL, Magnetic Resonance in Medicine 46:1-5 (2001), a magnetic 
field cycling method is described for transforming the proton spin-order to carbon spin 
polarization. In a first step the external magnetic field is reduced (from the relatively high 
geomagnetic field) bringing the combined proton-carbon spin system into its strong coupling 
regime. In this regime the scalar coupling (J-coupling) strongly influences the evolution of the 

30 spin system. The reduction of the field should be fast, giving a diabatic (non-adiabatic) 
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process. In a subsequent step the field strength is slowly increased (an adiabatic process). The 
field cycling will result in a substantial increase in the polarization of the spins of the 13 C 
nuclei, giving an in-phase NMR-signal, and increase the usefulness of the compound as a 
contrast agent for use in imaging procedures. However, the result of the field treatment will 
5 be dependent on the scalar coupling of the combined spin system and the properties of the 
magnetic field. In WO 00/71 166 examples of field cycling schemes are described that give a 
substantial improvement in the image quality. To make the method even more attractive for 
use in medical and diagnostic applications it would be of high value to further increase the 
degree of polarisation of the carbon spins and to prolong the relaxation time of the nuclear 
10 spin system. 

Summary of the Invention 

The object of the present invention is to. provide a method and an apparatus for producing 
MRI contrast agent with a high degree of polarisation of the imaging nuclei spins. 

The object is achieved by the method as defined in claim 1, the apparatus as defined in claim 
15 11, and by the computer program product as defined in claims 17 and 18. 

The method for producing contrast agent according to the present invention comprises the 
steps of obtaining a solution in a solvent of a hydrogenatable, unsaturated substrate compound 
and a catalyst for the hydrogenation of a substrate compound, hydrogenating the substrate 
with hydrogen gas (H2) enriched in para-hydrogen (p-'Hb) to form a hydrogenated contrast 
20 agent and exposing the contrast agent to a oscillating magnetic field in combination with a 
stationary magnetic field for enhancing the contrasting effects of the contrast agent adapted 
for use in an MR application. 

According to one embodiment of the present invention the oscillating magnetic field is 
oscillating with a frequency within the region of radio frequencies (around 10 Hz to several 
25 GHz) and preferably with a frequency in the interval 5 kHz to 500 MHz. 

According to a further embodiment of the present invention the step of exposure to the 
oscillating magnetic field in combination with the stationary magnetic field is performed 
during the step of hydrogenation, wherein the step of exposure is performed for reducing the 
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relaxation of the spin system of the contrast agent, whereby the contrasting effects of the 
contrast agent is enhanced. 

According to a still further embodiment of the invention the step of exposure to the oscillating 
magnetic field in combination with the stationary magnetic field is performed after the step of 
5 hydrogenation, the step of exposure is performed for enhancing the degree of polarization of 
imaging nuclei of the contrast agent, whereby the contrasting effects of the contrast agent is 
enhanced. The oscillating field is applied in the form of a series of pulses with the frequency 
and angle varying and adapted for the contrast agent in use. 

The apparatus for producing MR contrast agent according to the present invention comprises 
10 means for producing an oscillating magnetic field and means for producing a stationary 

magnetic field. The means for producing both the oscillating magnetic field and the stationary 
magnetic field may advantageously be of the same type as commonly used in NMR 
equipment. 

One advantage afforded by the apparatus and method according to the present invention is 
15 that contrast agent can be produced that significantly improves the image quality and/or 

speeds up the process in an MRI application and/or improves the analysis performance in an 
MNR application. 

A further advantage is that novel types of imaging, not possible, or very difficult to perform 
with prior art techniques, can be performed. 

20 The advantages are achieved by the apparatus and method according to the present invention 
by providing a high degree of polarization of the imaging nuclei of the contrast agent, by that 
the process of producing contrast agent is fast, thereby reducing the problems with relaxation 
of the spins system of the contrast agent and by that the relaxation of the spins system is 
effectively reduced. 

25 Brief Description of the Drawings 

The invention will now be described in detail with reference to the drawing figures, in which 



Fig. 1 is a flowchart illustrating the main steps of the method according to the invention; 
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Fig. 2 is a schematic drawing illustrating the apparatus according to the invention; 

Fig. 3 is a schematic drawing illustrating the magnetic treatment unit of the apparatus 
according to the invention is; 

Fig. 4 is a flowchart illustrating one embodiment of the method according to the invention; 
5 Fig 5 is a flowchart illustrating one embodiment of the method according to the invention; 

Detailed Description of the invention 

Parts of apparatus and parts of the process described in WO 00/71 166 are advantageously 
utilised also in the present invention. WO 00/71 166 teaches that the hydrogenation reaction is 
preferably performed by mixing gaseous para-hydrogen (or ortho-deuterium) enriched 
10 hydrogen with a solution of an unsaturated compound and a hydrogenation catalyst. 

The main steps of the of the method according to the invention will be described with 
reference to the flowchart of FIG. 1 and comprises the following main steps: 

100: obtaining a solution of a solvent, a hydrogenatable, unsaturated substrate compound and 
a catalyst for the hydrogenation of a substrate compound; 

15 110: introducing the solution into a chamber containing hydrogen gas (H 2 ) enriched in para- 
hydrogen (p- 'H2) in order to hydrogenate the substrate to form a hydrogenated contrast 
agent. 

The method according to the present invention introduces a main step of: 

120: subjecting the contrast agent to an oscillating magnetic field in combination with a 
20 stationary magnetic field. ; 

The oscillations of the oscillating magnetic field should preferably be in the kHz-MHz-region. 
Such oscillating fields are commonly known as rf-fields, and the term rf-field will hereafter 
be used to denote an oscillating magnetic field with a frequency in the interval 5 kHz to 500 
MHz. A frequency range of 10-500 kHz is preferred for many applications. The rf-field 
25 treatment enables spin-order to be transferred from protons in the freshly hydrogenated 
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contrast agent into polarization of a nucleus within the same molecule with a slower 
relaxation, preferably a ,3 C or ,5 N nuclei. These nuclei will be referred to as imaging nuclei. 
In addition an rf-field may be used to reduce the relaxation rate of the spin systems during the 
hydrogenation-process, giving a higher polarization of the contrast agent An apparatus and a 
5 method according to the present invention for providing the rf-field and stationary field 

during the production of the imaging agent is described below. The method and apparatus will 
be exemplified with contrast agent comprising carbon ( 13 C). This should be regarded as a non 
limiting example. Other isotopes could be used with slight modifications to the detailed steps 
of the method. 

10 The hydrogenatable substrate used may be a material such as a para-hydrogenation substrate 
as discussed in WO99/24080. For in vitro or in vivo MR studies of biological or quasi- 
biological processes or synthetic polymer (e.g. peptide, polynucleic acid etc.) syntheses, the 
substrate is preferably hydrogenatable to form a molecule participating in such reactions, e.g. 
an amino acid, a nucleic acid, a receptor-binding molecule, etc., either a natural such molecule 

15 or an analog. 

The solvent used in step 100 of the process of the invention may be any convenient material 
which serves as a solvent for the substrate and the hydrogenation catalyst. A number of 
possible solvents are discussed in WO99/24080. When the contrast agent is for use in in vivo 
MR investigations, the solvent is preferably physiologically tolerable. Water is a preferred 

20 choice of solvent, used in combination with a water-soluble catalyst. If other solvents that are 
not physiologically tolerable are used, the solvents has to be removed before use in a patient, 
for example by vacuum-spray. Other rapid solvent removal techniques, e.g. affinity 
techniques, may however be used. The solvent is preferably used at or near the minimum 
quantities required to maintain substrate, catalyst and contrast agent in solution without 

25 experiencing viscosity problems during the hydrogenation reaction. 

The hydrogenation catalyst is preferably a catalyst as discussed in WO99/24080, e.g. a metal 
complex, in particular a rhodium complex. 

The enriched hydrogen, which may be pure or 2 H2, or a mixture of ^ and 2 H2, may 
optionally contain other gases, but is preferably free from oxygen or other reactive or 



WO 2004/019995 PCTYEP2003/009587 



paramagnetic gases, may be prepared by cooling hydrogen, preferably to a temperature below 
80 K, more preferably to below 50 K, even more preferably to below 30 K and especially 
preferably to below 22 K, and allowing the nuclear spin states to equilibrate, optionally in the 
presence of a solid phase equilibration promoter, e.g. FejOij, Fe203, activated charcoal, etc. 
5 The enriched hydrogen is then preferably removed from the equilibrator and optionally stored 
before use. A method of preparation and storage of enriched hydrogen is described in 
WO99/24080. A preferred and novel method of storage and equipment for that purpose has 
been developed by the inventors. The enriched hydrogen is transferred to and stored in gas 
cylinders made of inert material. Inert in this context should be understood as made up by 

10 material essentially free from paramagnetic materials (primarily iron) and other para- 
hydrogen relaxing compounds (e.g. palladium). Examples of inert materials suitable for the 
gas cylinders are aluminum and carbon fiber reinforced epoxy. The enriched hydrogen will 
decay slowly (in the order of 10% per week). Such a decay rate is acceptable for most 
applications, especially compared with the cost and handling problems of previous storing 

15 methods involving storing at cryogenic temperatures. 

For the hydrogenation reaction, a reaction chamber is filled with enriched hydrogen 
optionally under pressure, preferably 5-20 bar, and the catalyst and substrate solution is 
introduced either as a thin jet, by spraying or by atomising, into this reactor. If desired, the 
solution may be produced by mixing separate solutions of catalyst and of substrate. To ensure 
20 proper mixing, a distributor or a plurality of spray nozzles may be used and the chamber 

contents may be mixed, e.g. by a mechanical stirrer or by appropriately shaping the chamber 
walls to promote turbulent mixing where there is a flow of reaction mixture in the chamber. 

The process may be performed continuously with a flow reactor, e.g. a loop or tube reactor, or 
alternatively it may be a batch-wise process. Preferably there will be a continuous or pulsed 

25 flow of enriched hydrogen and solution into the reactor, a continuous or batch-wise removal 
of liquid solution from the base of the reactor, and a continuous or batch-wise venting of 
unreacted gas from the reactor. The enriched hydrogen and solution passing into the reactor 
are preferably temperature-controlled to ensure the gas/droplet phase in the reactor is at the 
desired temperature. This can be achieved by providing input lines with temperature sensors 

30 and heating and/or cooling jackets. 
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If a non-aqueous solution has been used the contrast agent is preferably mixed with water 
after the hydrogenation and the exposure to the magnetic fields. The water used is preferably 
sterile and also preferably essentially free of paramagnetic contaminants. The resultant 
aqueous solution is then preferably treated to remove the hydrogenation catalyst, e.g. by 
5 passage through an ion exchange column, preferably one free of paramagnetic contaminants. 
The water may be temperature-controlled as may be a mixing chamber where water and 
contrast agent solutions are mixed so as to ensure the aqueous solution enters the ion 
exchange column at the appropriate temperature. Strongly acidic, sodium ion charged ion 
exchange resins such as DOWEX 1x2-400 (Dow Chemicals) and Amberlite IR-120 (both 
10 available from Aldrich Chemicals) resins may conveniently be used for the removal of typical 
metal complex hydrogenation catalysts. For fast ion exchange, the resin is preferably cross- 
linked to only a low degree, e.g. a 2% divinyl benzene cross-linked sulphonated, sodium ion 
loaded polystyrene resin. 

Removal of the non-aqueous solvent may then conveniently be effected by spray flash 
1 5 distillation e.g. by spraying the aqueous solution into a chamber, applying a vacuum, and 

driving the organic solvent free aqueous solution from the chamber using an inert, preferably 
non-paramagnetic gas, e.g. nitrogen. Indeed in general the flow of liquid components through 
the hydrogenation apparatus is preferably effected using applied nitrogen pressure, e.g. 2 to 10 
bar. 

20 The resulting aqueous contrast agent solution may be frozen and stored or may preferably be 
used directly in an MR imaging or spectroscopy procedure, optionally after dilution or 
addition of further solution components, e.g. pH modifiers, complexing agents, etc. Such 
direct use may for example involve continuous infusion or alternatively injection or infusion 
of one or more dose units. Bolus injection is particularly interesting. 

25 The whole process from beginning of hydrogenation to the delivery of the finished contrast 
agent in for example a syringe may conveniently be effected in less than 100 seconds, indeed 
it is feasible to produce dosage units in less than 10 seconds, which is substantially less than T\ 
for the potentially interesting imaging nuclei. 
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Preferably, the surfaces contacted by the contrast agent during the process of the invention are 
substantially free of paramagnetic materials, e.g. made of glasses as used for hyperpolarized 
^e containment as discussed in WO99/17304 or gold or polymer, optionally deuterated. 
Surfaces contacting a non-aqueous solvent (e.g. acetone) should be acetone-resistant and 
5 valves may be magnetically controlled and provided with solvent resistant Teflon or silicone 
parts. 

An apparatus suitable for producing contrast agent by the method according to the present 
invention will be described with reference to FIG. 2. Hydrogen (*H 2 ) enriched in para- 
hydrogen is fed from the para-hydrogen source 200 into a reactor 210. A hydrogenation 

10 catalyst solution from a catalyst reservoir 220 and a hydrogenatable substrate solution from a 
substrate reservoir 230 are fed into the reactor 210. The liquid settling in reactor 210 is 
transferred to a magnetic treatment unit 240, which essentially comprises a compartment for 
the solution surrounded by a low field, two-channel NMR spectrometer, and thence to 
finishing unit 250 for cleanup, quality control and possible addition of additives and solvent 

15 removal. The finishing unit may comprise an ion exchange column and a solvent removal 
chamber equipped with a spray nozzle. After the passage through the finishing units the 
imaging agent is delivered to, for example, a syringe for injection in a patient Alternatively 
the imaging agent is stored for later usage. The reactor 210, may in addition be provided with 
a initial magnetic treatment unit 260 adapted for producing a stationary magnetic field and an 

20 rf-field. As an alternative, the initial magnetic treatment unit 260 may be combined with the 
magnetic treatment unit 240 whereby all magnetic treatment will be performed within the 
reactor 210. 

An embodiment of a magnetic treatment unit 240 in accordance with the present invention 
will be described with reference to the schematic block diagram of FIG. 3, in which the main 

25 parts of the NMR-unit are outlined. The same type of NMR-unit may be utilized as the initial 
magnetic treatment unit 260, or a combined magnetic treatment unit. The NMR unit 
comprises a coil system 300, a radio frequency synthesizer 310, an amplifier 320, and a 
computer 330. The coils of the coil system 300 are connected to power supplies 345, which 
preferably are connected to and controllable by the computer 330. The coil system comprises 

30 a resistive magnet that gives a main field of 0.1-100 mT, preferably 1-10 mT. This magnet 
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may preferably be made as a Helmholz pair 340 with a large radius, around 30 cm, to give a 
high homogeneity over the sample space which preferably has a volume of 10-100 mL in the 
center of the coil system. The coil system 300 is preferably provided with three linear shim 
coils, x, y, and z (not shown) to allow for compensation of external field gradients, which 
5 primarily may be due to the presence of a nearby MRI system. The inhomogeneity of the 
Helmholtz pair will also be improved by the shims. Any other geometry (solenoidal) or 
magnet design (permanent magnet) compatible with the general design of the equipment can 
be used. Resistive magnets are suitable, due to the moderate field levels required. 
Superconducting magnets may also be used. Higher order shims can be added if there is need 

10 to place this equipment very close to the MRI magnet. When required, the magnet coils may, 
completely or in part, be wound from aluminum wire to reduce weight. Pairs of Hall plates 
(not shown) or other magnetic field detectors may also be present to allow for measurement 
of the field gradients in x, y and z directions so that the computer 330 can compensate for 
field disturbances automatically. The Hall sensors provide in addition the mean value of the 

15 magnetic field strength and is used to lock the field to the transmitter frequencies. 

An NMR-coil 360 is used to produce the excitation signal. The excitation signal is 
conveniently generated digitally. The rf-coil 360 is connected via switching means 365 and a 
amplifier 320 to a digital waveform generator 375, the amplifier and the waveform generator 
forming the rf-synthesizer 310. The digital waveform generator can be a stand-alone 
20 instrument or a PC card. The desired frequency spectrum and phase relations are created 
arbitrarily within the limits of the time resolution of the instrument/device. The excitation 
signal from the waveform generator is amplified to create the desired excitation field strength 
in the coil. Generally, the pulses will be long at the low frequencies anticipated, and 
consequently the excitation power is low, typically a few watts. 

25 The NMR coil 360 can be realized in a number of ways. The simplest is a solenoid, which has 
a high efficiency, and is a preferred geometry at low frequencies. The quality factor of the coil 
is maximized by the use of Litz wire. The coil may be tuned (double tuned) or non-tuned 
during excitation. Likewise, the coil can be tuned or non-timed in reception mode. The 
switching means 365 can be a mechanical relay or passive circuitry. 
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The receiver-part of the spectrometer is convenient for controlling and optimizing the 
spectrometer performance, but is not crucial for the production of the contrast agent. It can 
also be used for polarization measurement. An important part of the receiver-chain is the low 
noise amplifier 380, which is needed for enabling optimal performance of the spectrometer. 
5 The amplifier should have a sufficiently large bandwidth and low noise (e.g. less than 3dB). 
In the case that the coil is tuned the bandwidth may be increased by negative feedback (active 
damping). 

The amplified signal is digitalized directly and filtered to an appropriate bandwidth by an AID 
converter 385, preferably provided as a PC card. Interpolation and decimation reduce the 
10 signal frequency and number of sampling points. A Fast Fourier Transform (FFT) is applied 
to the signal in order to calculate the NMR spectrum. 

Preferably, most of the functions are integrated in a PC environment, in particular the 
frequency generation and sampling. The automatic shimming and field locking can also be 
handled from the PC. 

15 A method according to the invention, advantageously utilizing the above described apparatus, 
for increasing the polarization of a contrast agent, will now be described. The method causes 
a transfer of proton spin order to a J-coupled long-Ti nucleus by rf-pulse sequences. The 
pulses are applied to the contrast agent prior to injection into for example a patient The 
method is generally applicable for use with a large variety of substrates, but the pulse 

20 sequence will be specific for the substrate and the stationary magnetic field, Jfy. The method 
will be exemplified with a molecule containing only one none-zero spin nucleus, l3 C with 
S=l/2, which is hydrogenated with para-hydrogen and preferably in aqueous solution. It 
should be noted that the method according to the present invention is not limited to this 
example. The method of using rf-pulse sequences for spin order transfer is applicable to any 

25 molecules with spin 1/2 hetero nucleus. 



A prerequisite for the method is that the hydrogenation catalyst works by such a mechanism 
that the two hydrogen atoms from one molecule of para-hydrogen are delivered to one 
molecule of substrate to produce one molecule of product. These protons pairs are then 
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produced in a singlet spin state, which corresponds to a high degree of order, and the rf-pulse 
sequence is intended to transfer this order to carbon polarization. 

The hydrogenation is typically performed in the earth field, and the hydrogenation will occur 
at different times for the individual molecules. This will result in a system finally being in a 
5 statistical steady-state characterized by various populations of the energy states of the spin 
system, corresponding to the loss of the interstate initial correlations, and the subsequent loss 
of some of the initial order. The method according to the invention uses spin excitations by 
the way of radiofrequency pulses (rf-pulses) on either spin species, i.e. the frequency is 
adjusted to correspond to the Larmor-frequency of the nuclei. This is performed in a magnetic 

10 field sufficiently low such that the proton spins are in the so-called strong coupling regime, 
that is, the difference of the carbon-protons couplings and/or the difference of proton 
resonance frequencies due to different chemical shifts are not large compared with the inter- 
proton interactions. This has two important advantages. Firstly, there is no need to use a large 
magnetic field for the order-to-polarization step, and secondly the very existence of a strong 

15 inter-spin coupling is the root of the efficient manipulation of the spin state of the system. 

In the embodiments of the present invention the notations defining angles of the applied fields 
refers to a rotating frame and are the notations commonly used in the discipline of NMR. The 
phrases "pulse on carbon", "pulse on hydrogen" etc should be understood as a rf-pulse with 
the frequency adapted to the nucleus and the stationary field as to correspond to the Larmor 
20 frequency of the nucleus. 

Two alternative embodiments of the invention directly applicable in the case of the carbon 
nucleus having only two adjacent protons will be described. These prerequisites hold for e.g. 
2,3-dideuterosuccinate-l- ,3 C. The case when more than two protons are coupled to the carbon 
nuclues is more complicated than the former one giving a much more complex structure of 
25 the quantum energy states. This requires parameters such as the delay times between pulses 
and the angles of the applied fields to be adjusted by optimization, either by computer 
simulation or by experiments. 

The pulse sequence according to the first embodiment comprises two consecutive main parts: 
The first part is a preparation of the state of the system starting from the steady-state density 
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matrix at the end of the hydrogenation, through 180°-pulses followed by evolution periods, 
and the second part is production of a net spin S (carbon) polarization by way of non-180 0 
pulses (initially 90°) and evolution periods. In principle only the spins S are excited, but in 
practice one may apply 180° pulses on the spins / (protons) simultaneously with 180° pulses 
5 on the spins S to reduce the effects of field inhomogeneities. 

The steady-state density matrix, a u , after hydrogenation is given by the projection of the 
density matrix a = | (l - 4 i,i 2 ) on the eigen states of the Hamiltonian (rotating frame), 
H = J X X\* + JxJx,S, + > 

a * ~ i (* " 4 {fix 1 !* + P S i (A* J 7z ) + <AftJu + )))' The coefficients p and q depend on 
10 the scalar couplings, J ijt and determine the theoretical limit of the obtainable polarization. 

The first main part consists of transforming the steady-state density matrix above to a state 
that as closely as possibly resembles S x (l ly I 2x - I\ x h y )= S z K y y I to optimize the (scalar) 
constant in front of this operator. This is achieved by a series of 1 80° pulses and subsequent 
evolution periods. 

15 The second main part of the sequence starts with a 90y(5) pulse that converts S t K y into 
S x K y . The subsequent evolution transforms this state into 

2S x K y exp(- iHtp.S x K y exp(i7#) which is a linear combination of S y , 2K y S x and 
2K Z S X (where K x a \{l u -I 2x )). The coefficients in front of these terms oscillate, and the 
duration of the evolution will be chosen so that the term containing S y is maximized (at the 
20 same time 2^,^ vanishes). Then a 90 X (S) pulse is applied that produces a linear 

combination of S x and 2K y S x . Another (identical) evolution time gives a linear combination 
of the terms (S x + CS y ) and 2 K y S x . A pulse of angle ^ = arctan(C) , (f> x (5) , will give the 
terms S x and 2K y S x (the last term is unaffected by a pulse on S around Ox). Neglecting 
relaxation a subsequent evolution will not affect S x , whereas the term 2K y S x transforms as 
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above. Repeating the above procedure, changing the angle 0 accordingly, increases the 
polarization of carbon (i.e. S z ). 

A realization of the first embodiment of the method according to the invention will be 
described with reference to the flowchart of FIG. 4. The realization of the first embodiment 
comprises the steps of: 

400: Applying a stationary magnetic field Bo. 

410: Hydrogenation of the substrate with para-hydrogen in the presence of the magnetic 
field B 0 . 

420: Applying a series of 1 80° x pulses followed by delays t§ on carbon. The delays after 
each pulse will typically differ. The delays are typically in the order of 1-100 ms. 
The number of repetitions are typically one or two but may be more. The purpose of 
this step is to bring the system into a state consisting of a zero quantum coherence 
involving the two protons and the carbon. The determination of delays and number of 
repetitions can be made analytically (if having one carbon and two protons) or by 
numerical methods (more spins) and involves maximizing the term S^Ij^urhx^y in 
which S refers to the spin operator of carbon and Ij and I 2 refer to the spin operators of 
the two protons. Step 420 corresponds to the first main part of the embodiment 

430: Applying a 90° y pulse on carbon. 

440: Waiting for til s. 

450: Optionally applying simultaneous 1 80° x pulses on hydrogen and carbon in order to 
compensate for the effect of field inhomogeneities (spin echo pulses). 

460: Optionally waiting for t/2 s. 

470: Applying a pulse fa on carbon, the angle fa may be different in each step and as 

previously described determined by analytical calculation or by numerical simulation. 
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480: Repeating steps 440 to 470 m times. Steps 440-470 gives a progressive build up of 
carbon polarization in the direction of the external field axis. The number of 
repetitions, m 9 is typically 1-5. The steps 430-480 corresponds to the second main part 
of the embodiment. 

5 The first embodiment is characterized in that the pulses substantially effecting the polarization 
are applied only to the carbon nucleus, i.e. with a frequency corresponding to the Larmor 
frequency of carbon. The step involving pulses on both carbon and hydrogen is only taken to 
compensate for experimental imperfections. Depending on the molecule, a carbon spin 
polarization between 50% and 85% is typically achieved. The total duration of the 
10 polarization process is about one to several hundred milliseconds. Thus, the method is fast 
which can be crucial because of the loss of polarization due to spin relaxation. 

The second embodiment comprises three consecutive main parts. The first main part 
transforms the steady-state density matrix (see above) to a state that as closely as possibly 
resembles {(l -4{l u I 22 -I u I 2x -^ 2y ))- ifr-^C-Va* i e - optimizes the (scalar) 

15 constant in front of the operator K x (optimally equal to one). This is accomplished by a series 

of 180° pulses and subsequent evolution periods of equal duration. At the end of this part the 
density matrix is proportional to (disregarding small terms and the unimportant constant term) 
I u I 2x + cK x (c being close to unity). The second main part starts with a 90y(/) pulse that 
produces the term [l u I 2x - Iiyliy)- A subsequent evolution period gives the term 
20 S x {l ly I 2x + Juhy )• The main with a 90y(5) pulse transforming S x -» S x and a 

90^ (i) pulse, with a suitable phase, q> , that transforms {l X7 I 2x + A* ) * nto ^2x ~ ^\ y h y )• 
The useful part of the density operator is now 2SJ lt I 22 . An evolution period produces the 
useful part S x (neglecting terms that are of no interest) which is finally transformed to S x by 
a 90y(S) pulse. 

25 A realization of the second embodiment of the method according to the invention will be 

described with reference to the flowchart of FIG. 5. The realization of the second embodiment 
comprises the steps of: 
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500: Applying a stationary magnetic field Bo- 

510: Hydrogenation of the substrate with para-hydrogen in the presence of the magnetic 
field B 0 . 

520: Applying a series of 1 80° x pulses followed by delays // on caibon. The number of 
repetitions are typically one or two but may be more. The purpose of this step is to 
bring the system into a state consisting of a two proton zero quantum coherence. The 
determination of delays and number of repetitions can be made analytically (one 
carbon, two protons) or by numerical methods (more spins) and involves maximizing 
the factor (IiJ2x-Iiyl2y). Step 520 corresponds to the first main part of the second 
embodiment. 

530: Applying a 90° y pulse on hydrogen, which results in a two-proton-double quantum 
coherence. 

540: Waiting for fc/2 s. During this time the two-proton-double quantum coherence 

transforms into a three-spin coherence involving the carbon spin. Steps 530 and 540 
corresponds to the second main part of the second embodiment. 

550: Optionally applying simultaneous 180° x pulses on hydrogen and carbon in order to 
compensate for the effect of field inhomogeneities (spin echo pulses). 

560: Optionally waiting for t 2 l2 s. 

570: Applying simultaneous 90° y and 90% pulses on carbon and hydrogen, respectively, for 
producing a transverse carbon polarization. 

580: Waiting for 0/2 s. 

585: Applying simultaneous 1 80° x pulses on carbon and hydrogen. 
590: Waiting for /j/2 s. 
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595: Applying a -90° y pulse on carbon, for transforming the transverse carbon polarization 
into a longitudinal polarization, which is parallel to the applied field. Steps 570-595 
corresponds to the third main part of the second embodiment. 

It should be noted that in the second embodiment, contrary to the first embodiment, rf-pulses 
directly effecting the polarization are applied to both the carbon and hydrogen nuclei. The 
polarizations achievable by the second embodiment are typically between 60% and 98%, 
depending on the molecule. For most instances, they are larger than in the first embodiment, 
but the pulse sequence is in most cases more complex and the duration of the polarization 
sequence is in general longer. Hence, the loss due to spin relaxation during the increased 
process duration has to be considered. Which of the embodiments to use under practical 
conditions is thus decided based on these factors and on the capability of the equipment. The 
decision is typically "case by case". Implementation examples and simulated results are 
enclosed in Appendix A. As appreciated by the one skilled in NMR and MRI, the tools and 
practices utilized in determining the pulse sequences in NMR and MRI applications can also 
be used as a guidance for determination of the pulse sequences used for producing contrast 
agent according to the present invention. 

It is conceivable to use pulses that are modulated in amplitude, phase or frequency or any 
combination of those. Alternatively, the individual pulses can be substituted with composite 
pulses (i.e. sequences of phase-shifted rf pulses). This provides better phase-control over a 
20 wider bandwidth. This gives better results in practice and relaxes the demands on precision of 
the hardware. Examples of composite pulses can be found in: "Band-Selective 
Radiofrequency Pulses", Geen, H. Freeman, R. J. Magn. Res. 93, 93-141(1991), "Composite 
Pulses" Levitt, M. Progress in NMR spectroscopy* 18, 61-122 (1986), "Composite Pulses 
without Phase Distortion" Tycko, R. J. Magn. Res. 61, 90-101 (1985) and "Broadband, 
25 Narrowband, and Passband Composite Pulses for Use in Advanced NMR Experiments" 
Wimperis, S. J. Magn. Res. SerA 109, 221-231 (1994). 

In both the described embodiments the delay times, //, are typically in the order of 1-100 ms. 
The amplitude of the rf-fields is typically 1-100 /xT, the duration of the pulses follows from 
the rf-amplitude and the frequency adapted to the Larmor frequency of the target nucleus. The 
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stationary magnetic field B 0 , applied in step 400 and 500 respectively, is preferably 0.1-100 
mT, and more preferably 1-10 mT. 

The method according to the invention will facilitate the production of highly polarized 
contrast agent. The properties of the contrast agent may further be improved if the relaxation 
5 rate of the spin system is reduced. One way of reducing the relaxation rate of the spin system, 
i.e. extending the useful lifetime of the final contrast agent, is by decoupling of the former 
para-hydrogen protons during the process of hydrogenation. The apparatus and the method 
according to the invention of subjecting the contrast agent to rf-pulses, may advantageously 
be used also to reduce the relaxation of the spin system. 

10 During the process of hydrogenation a substantial relaxation of the spin system occurs. The 
relaxation takes place both when the para-hydrogen is attached to the catalyst in the 
intermediate state and later, when it is transferred to the molecules of the contrast agent. The 
relaxation results in a reduced PHTP effect and impairs the performance of the contrast agent. 
During the hydrogenation different molecules will be hydrogenated at different times and the 

15 fonner para-hydrogen protons will no longer be in the singlet state. Intramolecular dipole- 
dipole coupling will cause relaxation of the spin system. The dipole-dipole interaction 
contribution to the relaxation can be significantly reduced by forcing the protons to remain in 
the singlet state. From a quantum mechanical consideration of the spin system it can be 
realized that, by letting the process of hydrogenation take place in the presence of both a rf- 

20 field and a static magnetic field, the protons can be forced to remain in their singlet state. 
Hence the relaxation can be significantly reduced. 

In order to utilize the method of exposure to rf-field and magnetic field during hydrogenation 
the reactor 210 has to be equipped with means for producing the fields. In addition the 
material of the reactor 210 has to be transparent both to the rf irradiation and the static 

25 magnetic field. Preferably a non-conducting, i.e non-metallic material is used, for example 
glass-fibre reinforced epoxy. This material is also capable of withstanding the substantial 
pressure used during the process. The external static magnetic field and the rf-fields can 
preferably be provided by a magnet system similar to that of the magnetic treatment unit 240 
described above. Alternatively the external magnetic field may be provided by a solenoid coil. 

30 The required field strength is in the order of 1 mT and the requirement of homogeneity is 



WO 2004/019995 



19 



PCT/EP2003/009587 



fairly moderate - inhomogeneities of up to ±6% are acceptable. The rf-field may preferably be 
provided by crossed orthogonal saddle coils, which can produce circularly polarized rf-fields 
as well as more complex decoupling sequences. The frequency of the rf-field should 
correspond to the Larmor frequency of the protons, e.g. 42.6 kHz if the external field is 1 mT 

( v = > where vis frequency of the rf-field, y the gyromagnetic ratio and B 0 the magnetic 

field). The intensity of the rf-field should preferably be optimized for each case. A suitable rf- 
field amplitude will typically be approximately 50 /*T, with an acceptable inhomogeneity of 
the order of 10 pT. The generation and control of the required fields are well-known in the 
area of NMR and the chosen implementation can be modified in many ways as appreciated by 
those skilled in the art. All the magnetic field treatments are preferably done with the same 
magnet system and in the same chamber. Thus, in this preferred embodiment of the apparatus, 
the reactor 210 and the magnetic treatment unit 240 are combined into one unit. 

The embodiments of the method according to the invention described above may be modified 
to also include the magnetic treatment for reducing the relaxation by letting the step of 
hydrogenation (step 400, 500) to be performed in the presence of a combined external 
magnetic field and an rf-field. The step in this preferred embodiment being: 

400, 500: Hydrogenation of the substrate with para-hydrogen, performed in the presence 
of an external stationary magnetic field and an external rf-field. 

The method of exposing the contrast agent to a combined external magnetic field and an rf- 
field during the hydrogenation may be combined also with methods of spin order transfer 
other than the here described method utilizing series of rf-pulses. For example, the method of 
using field cycling schemes described in WO 00/71 166, may be improved if a combined 
external magnetic field and an rf-field is present during the hydrogenation. 

As the magnetic treatment unit 240 in a preferred embodiment of the invention has all the 
facilities of a two-channel NMR spectrometer, it can be advantageously used for 
characterizing and checking the quality of the produced contrast agent. The ability of 
characterizing, for example measuring the degree of polarization in the produced contrast 
agent, is valuable in developing and fine-tuning the pulse sequences of the above-described 
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embodiments. This is obviously most useful in the cases when the parameters of the pulse 
sequences have to be established from experiments, but also if the parameters are determined 
analytically or by computerized optimization, the degree of polarization may be increased by 
fine-tuning parameters, for example the delay times, experimentally. 

5 The method of the invention is preferably implemented by means of a computer program 
product comprising the software code means for performing the steps of the method by 
controlling parts of the apparatus according to the invention. The computer program product 
is typically executed on the computer 330. The computer program is loaded directly or from a 
computer usable medium, such as a floppy disc, a CD, the Internet etc. 

1 0 From the invention thus described, it will be obvious that the invention may be varied in 

many ways. Such variations are not to be regarded as a departure from the inventive concept, 
and all such modifications as would be obvious to one skilled in the art are intended for 
inclusion within the scope of the following claims. 
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Appendix A 

Exemplary pulse sequences: 

Pulses given in degrees with phase as subscript, delays as time t with index as subscript. 
Sequence 1-2 corresponds to the first embodiment of the invention, sequence 3 to the second 
5 embodiment and sequence 4 is a special case related to double hydrogenation. 

Sequence la: 

13 C 180 x -t,-180 x - t 2 -90 y -t3/2-180 x -t3/2- <f>\ x 
«H 180 x 

Sequence lb: 

10 l3 C 180 x -t,-180 x -t2-90 y -t3/2-180 x .t3/2- 01 R - 1^2-180*- ta/2.- <£2 X 
*H 1 80 x -t 3 1 80, 

Sequence lc: 

l3 C 180,-t r 180,- t 2 -90 y -t3^-180 x -t3^- <f>\ % - t3/2-180 x - \Jl r 02,- tjl -180, - tjl- 03 x 
l H . 180,-t 3 180,-ta 180, 

15 Sequence Id: 

,3 C 180 x -t,-180 x - t r 90 y -ty2-l 8(^-13^- <f>l x - t3/2-180 x - X^2 r 02 x - XJ2 -180,- XJ2- 03,- t3/2-180 x - tj/2- 04 x 
»H 180 x -t3 180 x -t 3 180,-t 3 180, 

Sequence 2a: 

*C 180,-t 3 -180 x -t,-180 x - 12-90^3/2-180^5/2- if>\ x 
20 'H 180, 

Sequence 2b: 

,3 C 180 x -t 3 -180 x -t r 180 x - 12-90^-180^3/2- 01,- Xjl-lZO* - Xyf2 r 02 x 
»H 180,-t 3 180 x 

Sequence 2c: 

25 l3 C 180,-t 3 -180 x -t,-180 x - 12-90^13^-180,^3/2- 01,- tj/2-180, - tj/2- 02,- xJ2 -180, - ta/2- 03 x 
l H 1 80,-13 1 80,-t 3 180, 
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Sequence 3: 

l3 C (180,-t,)^ 180^2/2-90^3/2-180.- tjtf^Oy 

'H 90 y -t 2 /2- 180^2/2-90^ -13/2-180* 

Sequence 4: 

,3 C 180,-1,-90^2^.180^12/2-90, 
'H 180 x 



Examples of predictions 

Contrast agent: Maleic acid 
10 Jhh=10.65 Hz, J H ic=15.5 Hz, Jh2C=05 Hz 

Pulse Sequence la, B 0 =5 mT, t,= 20.76 ms, t 2 =33.91 ms, t 3 =38.39 ms, 01=90°: 
Polarization P= 76.98 % , Total duration = 93.06 ms. 

1 5 Pulse Sequence lb, B 0 =5 mT, t,= 20.76 ms, t 2 =33.91 ms, t 3 =38.39 ms, 01=90°, 02=- 1 8.62°: 
Polarization P= -81.23% , Total duration = 131.45 ms. 

Pulse Sequence lc, Bo=5 mT, t t = 20.76 ms, t 2 =33.91 ms, t 3 =38.39 ms, 01=90°, 02=-18.62°, 
03=6.14°: 

20 Polarization P= 81 .7% , Total duration = 169.83 ms. 

Pulse Sequence 3, Bo=5 mT, n=2, ti=38.4 ms, t 2 =31.25 ms; t 3 =69 ms: 
Polarization P= 85.4% , Total duration = 177.0 ms. 



25 Contrast agent: succinic acid (2.3-D?) 

J H Hvic=7 Hz, J H ic=7 Hz, Jh2C= 0 Hz, A5H= 0 ppm 

Pulse Sequence la, B 0 =5 mT, ti= 44.09 ms, t 2 =46.64 ms, t 3 =63.89 ms, 01=90°: 
Polarization P= 71.55 % , Total duration = 154.62 ms. 
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Pulse Sequence lb, B 0 =5 mT, t,= 44.09 ms, t2=46.64 ms, t 3 =63.89 ms, 01=90°, 02=-3O.96°: 
Polarization P= -83.45% , Total duration = 21 8.5 1 ms. 

Pulse Sequence 1c, B<f=5 mT, tj= 44.09 ms, t 2 =46.64 ms, t 3 =63.89 ms, 01=90°, 02=-3O.96°, 
03=17.16°: 

Polarization P= 87.33% , Total duration = 282.4 ms. 

Pulse Sequence Id, B 0 =5 mT, ti= 44.09 ms, t 2 =46.64 ms, t 3 =63.89 ms, 01=90°, 02=-3O.96°, 
03=17.16° 03= -10.04°: 

Polarization P= -88.69% , Total duration = 346.29 ms. 

Pulse Sequence 3, Bo=5 mT, n=3, ti=63.9 ms, t 2 =71.4 ms, t3= 135.4 ms: 
Polarization P= 93.2% , Total duration = 398 ms. 
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Contrast agent: OAcetvl lactic acid (3,3-D?) 
J HH =7.1 Hz, J H ic=7.3 Hz, Jh2C=2.4 Hz, A5H= 3.2 ppm 

Pulse Sequence 2a, Bo=S mT, ti= 40.07 ms, t 2 =58.12 ms, t 3 =66.57 ms, 01=90°: 
5 Polarization P= 57.6 % , Total duration = 23 1 .33 ms. 

Pulse Sequence 2b, B<f=5 mT, t]= 40.07 ms, t 2 =58.12 ms, t 3 =66.57 ms, 01=90°, 02=-38.2°: 
Polarization P= -72.06% , Total duration = 297.9 ms. 

10 Pulse Sequence 2c, B 0 =5 mT, ti= 40.07 ms, t 2 =58.12 ms, t 3 =66.57 ms, 01=90°, 02=-38.2°, 
03=25.96°: 

Polarization P= 78.80% , Total duration = 364.47 ms. 

Pulse Sequence 3, Bo=5 mT, n=4, t|=63.9 ms, t2=71.4 ms, t3=135.4 ms: 
1 5 Polarization P= 93.07% , Total duration = 424.7 ms. 

Special case:Succinic acid, produced from acetylene dicarfaoxvlic acid bv two parahvdrogen 
reactions in rapid sequence: 
20 JHHvic=7 Hz, JHH g em=14 Hz, J H ic=7 Hz, Jh2c= 0 Hz, A8H= 0 ppm 



Pulse sequence 4, ti=28.43 ms, t 2 =56.9 ms: 
Polarization P=57% , Total duration =85.3 ms. 



